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^"^ ' Abstract. We study the experimental signatures of several states of a Bose-Einstein 

condensate of spin-1 atoms by quantum trajectory simulations of Stern-Gerlach 

^^-H I experiments. The measurement process itself creates an apparent random alignment 

jy-j ' for the spins, so that it is difhcult to distinguish between a condensate that initially has 

an alignment in an unknown direction and one with no alignment at all. Nonetheless, 
in repeated experiments with identically prepared condensates, it is possible to 
discriminate between certain ground states of the condensate proposed in the literature. 
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1. Introduction 

With tlie introduction of the optical trap in the experiments on Bose-Einstein 
condensation, simuhaneous confinement of atoms in several angular momentum states 
is now feasible |I[]. The atomic spins interact in spin-exchange collisions, which gives rise 
to collective many-spin states. Energy differences between several conceivable ground 
states of the spinor condensate p|, |, ^ have turned out to be small, so that it need not 
be obvious which one gets prepared in an experiment. 

On the other hand, resolving the contradicting opinions of the past P, ^, recent 
analyses have shown 0, §, |^ that a Bose-Einstein condensate may display signatures 
of spontaneously broken gauge symmetry, such as interference of two independently 



prepared condensates [jTO| or the Josephson effect, even if the gauge symmetry of the 
state is in fact not broken. The lesson is that the relation between the state of the 
system and the measurement results may be conditioned as much by the experiment as 
by the state itself. How to distinguish experimentally between the proposed states of 
the spinor condensate 0, |] is therefore a nontrivial question in its own right. 

In this paper we study the experimental signatures of a spin-1 condensate for three 
different initial states, "coherent state" with spontaneous alignment 0, "fragmented 
state" in which there are two macroscopically occupied spin states 0], and rotationally 
invariant "singlet state" 0]. Our main tool is quantum trajectory simulations [jTT], H, ^ 
of repeated Stern-Gerlach experiments. Once more 0, ^, 0], the interplay between 
state and measurement is evident. It appears that the coherent and fragmented states 
are experimentally indistinguishable, but can be distinguished easily from the singlet 
state. We also present a testable prediction for the distribution of the outcomes of the 
Stern-Gerlach experiments that probes the statistic of spin projections to all orders, 
not just averages and standard deviations. Finally, as an important by-product of our 
development we note that, while the simulations at first sight only apply to a quite 
artificial experimental scheme, they implement what we believe is a universal theory for 
spin measurements on a condensate. 

2. States of spinor condensate 

We consider putative ground states for the spinor condensate in the case when the 
spin species do not separate, but the atoms with the three angular momentum 
projections m = ±1 and m = commingle with identical spatial wave functions. 
We assume scattering lengths such that the interaction between the spins is effectively 
antiferromagnetic. For brevity we assume zero magnetic field. Nonetheless, we still pick 
a quantization axis z. We write the annihilation operators for the atoms with the three 
spin components along the quantization axis as a± and Qq. For notational simplicity, 
we take the total number of atoms A^ to be even. 

The coherent state is arrived at via the semiclassical argument, which treats the 
annihilation operators a± and qq as c numbers Q. The c numbers come with phases 
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that are not uniquely determined by minimization of energy. Assigning a value for the 
phases is tantamount to inserting by hand a form of spontaneous symmetry breaking. 
The corresponding quantized ansatz for the state of the spinor condensate reads |^ 

\C)= ,^ =(e-^V +e'%U^|vac). (1) 

Positing the state \C) is analogous to assuming that a single-component condensate is 
in a coherent state, as opposed to, say, a number state of the atoms 0, |], ||. x is an 
angle that characterizes the symmetry breaking. It selects preferred directions in the 
xy plane. 

The fragmented or "coherent-fragmented state" [Q is just a number state with half 
of the spins in the state + and half in the state — , 

|F) = |iV/2,0,iV/2). (2) 

The arguments give the numbers of atoms with the z components of the angular 
momentum equal to —1, 0, and +1. In the thermodynamic limit the state \F) should 
model the ground state in the presence of even the most minute magnetic field in the z 
direction |Q. While this state could behave differently in the x and z directions, it has 
no built-in structure in the xy plane. 

The singlet state of is of the form 

N/2 

\S) = Y.Ak\k,N-2k,k). (3) 

fc=0 

The coefficients A^. satisfy the recursion relation 



N-2k + 2 ^ .^. 

^^ = -Viv^2^TI^^-^' ^'^ 

and are then fixed by normalization, except for an arbitrary overall phase. The state 
\S) is the unique ground state for exactly zero magnetic field. It has the same form in 
all rotated frames, and thus possesses no intrinsic direction. 

3. Measurements on spins 

In our model for the measurements, we assume for the time being that one atom at a 

time is removed from the condensate and is subject to a Stern-Gerlach experiment that 

probes the spin component in the xy plane in a fixed direction at an angle ip with respect 

to the X axis. The observed sequence of spin projections is recorded and analyzed. All 

measurements are supposed to be completed in a time shorter than the evolution time 

scale of the spins in the condensate owing to spin-spin interactions. 

From angular momentum algebra it is easy to see that the annihilation operators 

for the spin states in the direction ip are given in terms of the original quantization axis 

operators as 

gj'/' 1 e"**^ 

a±{^) = —a^±-^ao + ^-a+, (5) 
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We carry out conventional quantum trajectory simulations of the measurement 
sequence [H, || numerically. Thus, suppose that, entering the nth measurement, the 
normalized state of the spins is \ipn)- We first calculate the probabilities for each spin 
projection, 

pin) ^ {'^n\aL{v)am{v)\i^n) ,„s 

N-n ^^ 

for m = — 1, and +1. Second, we use a random number generator to pick one of the m 
values in such a way that the probability for choosing m equals -P^-*. Third, we reduce 
the wave packet according to the selected value of m, so that the state vector for the 
(n + l)th measurement is 

i^„..)^ , °';'^"'""> . (8) 

The choices of the spin projection m for each measurement n constitute the data. 

We first consider the coherent state, \ipo) = \C). It turns out that, except for 
an inconsequential overall phase, and independently of the past choices of the spin 
projections, the state vector {ipn) is always a coherent state of the same form as \C), 
except that of course the total number of atoms is N — n not N. Furthermore, the 
probabilities for the outcomes of the Stern-Gerlach experiment are always the same, 

Pi") = icos2(v.-x), Po^"^ = sin2(^-x)- (9) 

Successive Stern-Gerlach experiments are therefore uncorrelated. The nature of the 
symmetry breaking is alignment. Namely, in the measurement directions (p = x^\t^ 
one only sees the result m = 0, whereas in the directions ^ = x ^"^^ V^ = X + ^ oii6 finds 
m = ±1 with equal probabilities, and no m = 0. 

Now, according to the notion of spontaneous symmetry breaking, the angle x varies 
at random from one spinor condensate to the next, and the experimenter has no a priori 
way of knowing it. As the observed frequencies of spin projections Nq, N± vary wildly 
with the angle x? measurements on just a single condensate do not seem to make a 
particularly discriminating test for or against the state \C). 

To gather more incisive data, we imagine repeating the experiment with a large 
number of identically prepared condensates. In principle we should discuss the 
frequencies of the observed spin projections, but in our simulations we have access 
to the probabilities Pq and PJ" as well. For brevity, here we often pick a particular 
measurement n, and focus on the combination of the probabilities 

p = pin) + pM _ pin) (^0) 

For a given condensate in the coherent state \C) this has the value P = cos'^{ip — x) ^ 
sin^((y9 — x), but the angle x is unknown. Nonetheless, if the angle x is evenly distributed 
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over the unit circle [0, 27r), we may calculate the probability density f{P) for the values 
of P. The result is 



1 



^e(-l,l); 



/(P) = <^ 7rVl-P2' ' ^' (11) 

[ 0, otherwise . 

Thus, we envisage using a large number of condensates. For each condensate one 
extracts the combination of probabilities P, in our numerical experiments using (p!o|) 
for some given measurement n, and in real experiments by estimating the probabilities 
using the observed frequencies of spin projections for each condensate. Finally one 
compares the distribution of the numbers P with predictions such as (|Tl|) . In the case 
of the coherent state \C), the numerical experiments simply amount to testing of our 
algorithms. 

We next move on to the fragmented state \F). In figure |^(a) we plot the numerically 
simulated probabilities Pq and PJ" as a function of the measurement number n for 
a sample condensate with A^ = 1000 atoms, and figure |I](b) gives the corresponding 
cumulative frequencies A^o^ N±- Notably, after an initial transient, the probabilities 
approximately stabilize at values compatible with the observed relative frequencies of 
the spin components. 

The expectation values for the numbers of spins with given projections are 
{al.{ip)a±{(p)) = ^N, {al{(p)ao{ip)) = ^N. One might naively expect that, even in a 
single condensate, one would find the spins approximately in the ratios N_ : Nq : iV_|_ = 
1:2:1. This clearly is not the case. In figure |I|(b) the ratios rather are 2:1:2. Moreover, 
the observed ratios vary at random from one condensate to another. 

The measurements are correlated; the outcome of an observation of the projection 
of a spin affects the state, which in turns affects the outcome of future observations. The 
correlations work out in an interesting way. In the first measurement the probabilities 
are \ '■ ^ '■ \, but then they quickly drift away and stabilize at some other values. Once 
the probabilities have stabilized, Stern-Gerlach measurements seem to be approximately 
independent repetitions of one another. However, the apparent state that is being 
repeatedly measured varies at random from one condensate to the other. Measurements 
on the fragmented state \F) in an individual condensate behave remarkably like 
measurements on the symmetry-broken coherent state \C). 

It is legitimate to ask if the states \F) and \C) can be distinguished at all. We 
address this question further in figure |^. We take 1000 spinor condensates, each with 
A^ = 500 atoms. For each condensate we record the quantity P of ( p!OD for the Stern- 
Gerlach measurement number n = 101. We then bin the results into 40 equally wide 
slots of P, and draw the histogram as circles. Also shown as the solid line is the 
prediction for the histogram if the probability distribution for the values of P is given 
by ([Tl|), as appropriate for the coherent state \C). Given the statistical fluctuations, we 
cannot tell the states \F) and \C) from one another. 

The situation is remarkably similar to the one that was encountered earlier with a 
single-component condensate when it came to the question of a number state, versus a 
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Figure 1. Simulated probabilities (a) and cumulative frequencies (b) for spin 
projections as a function of measurement number in one particular condensate that 
starts out with 1000 atoms in the fragmented state \F). 
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Figure 2. Histogram of the values P as in (Hfl) from numerical simulations of the 
fragmented state \F) (circles), and as deduced from (ni]) for the coherent state \C) 
(solid line). 1000 independently prepared condensates with N = 500 atoms each are 
used in the simulations. The range (—1, 1) of the quantity P is divided into 40 equally 
wide bins, and the vertical axis shows the number of P values in each bin. We use the 
measurement n = 101 in the simulation data. 



coherent state whose phase arises from spontaneously broken gauge symmetry [0, §, ^. 
When detection of the atoms is exphcitly considered, it turns out that both states lead 
to interference between independently prepared condensates, even though the number 
state seemingly does not provide any phase for the interference 0, ^. Both coherent- 
state and number-state condensates even exhibit the Josephson effect [^, Q. It is 
possible in principle to tell the difference between a coherent state and a number state 



by studying atom statistics for a small number (at most tens) of atoms |T^, but an 
experiment discriminating between coherent and number states is yet to be carried out 
in a condensate of any size. We believe that, analogously, there would be little or no 
practical difference between the states \C) and \F) of the spinor condensate. 

Let us finally consider the singlet state \S). The expectation values of the spin 
components are all {al.{(p)a±{(p)) = {al{(p)aQ{(p)) = ^N. However, just like in the case 
of the fragmented state \F), in a single condensate both the simulated probabilities and 
the frequencies stabilize at values that may be completely different from A^_ : A'^o : A^+ 
= 1:1:1. One may, of course, draw a histogram such as in figure ^ also for the state 
I S) , which we have done in figure ^ The discrepancy with the prediction of the coherent 
state \C) is evident even without any statistical analysis. 

Drawing the histogram for the state IS") might already have been be a case of 
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Figure 3. Same as figure eI except that the singlet state \S) was used in the 
simulations. 



overkill. While the frequencies for spin projections are seemingly random for a single 
condensate, we have found that they behave as expected when averaged over many 
condensates: 1 : 2 : 1 for the states \C) and \F), 1 : 1 : 1 for the state 15). The 
singlet IS*) should be experimentally distinguishable from the other two states simply 
by averaging the frequencies of spin projections over many condensates. 



4. Comparing schemes for Stern-Gerlach experiments 

Up to now we have envisaged picking one atom at a time out of the condensate. In 
reality one might rather do an en masse Stern-Gerlach experiment, subjecting the entire 
spinor condensate at once to a gradient of the magnetic field. The spin projections then 
separate in space, and can be counted by counting atoms. 

We assume that the detection is in principle still one atom at a time, as if the atoms 
fell on small detectors like photons hitting silver bromide grains on a photographic film. 
It is known from the theory of photon detectors |T^ that the joint counting rate of an 
array of absorbing broadband detectors is proportional to normally ordered correlation 
functions of electric field operators. Analogously, we conjecture that the joint probability 
of seeing a sequence of spin projections mi,m2, ■ ■ ■ , ruN is given by 

P(mi, . . . ,m7v) = K{a^^^{^) . . . al^^{(f)a^^{(f) . ..amM) , (12) 

where K is the normalization constant needed to make the sum of all probabilities 
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equal to unity. In photon detection theory the electric field operators are in time order 
as well, but inasmuch as we assume that we may ignore the reversible evolution of the 
spins during the measurement, the annihilation (creation) operators in (|12]) are taken 
at the same time and commute. The joint probabilities are then independent of the 
order of the indices mi, . . . ,mN, i-e., independent of the order in which the atoms are 
detected. 

Suppose that the spin system starts out with the state vector {ipo), and consider 
the correlation functions as in (|l2l). We may, of course, write 

(V'o|a]ni(<^) • ..al^{^)am^{^) . . . amAv)\ipo) 

Trivial as this rewrite is, it immediately leads to two crucial observations. First, in (P^) 
the probabilities for the various outcomes mi are precisely as given by (|^ for n = 1. 
Second, once any particular value mi has been picked, the process starts over; we are 
left with the correlation functions 

ii^l I «m2 (<^) • • • «mjv (<^)«rnjv (<^) • • • (^m^ (<^) \i^l) 

where the new state \ipi) is the n = 1 version of (P). We have thus shown that 
our quantum trajectory simulations produce spin sequences mi, . . . ,mN that have 
the probabihties predicted by (|T^). A fortiori, the results from quantum trajectory 
simulations agree with the predictions of (|1^). 

Our simulations were seemingly for a scheme in which one atom at a time is 
subject to a Stern-Gerlach experiment. We did not offer any concrete suggestions 
on how to carry out such experiments, but neither did we have to. Inasmuch as our 
reasoning leading to (|1^) is correct, our simulations also model the usual |]1| Stern- 
Gerlach measurements on the spinor condensate. 

One might also study the number operators for different spin projections, such 
as Uoif) = al{ip)ao{^), in total disregard of the measurement process. As far as we 
can tell, our quantum trajectory simulations averaged over many condensates give the 
same qualitative and semi-quantitative results that one would expect on the basis of 
elementary quantum mechanics with the operators uoif), n±{(f). 

As a qualitative example, the fluctuations oino{ip) in the state \F) are large, Ahq ~ 
a/| N. In quantum trajectory simulations one correspondingly sees wildly varying 
ratios in the frequencies of spin projections between individual condensates. From this 
viewpoint, measurement-induced alignment is the reason for the large fluctuations in 
no{(p) [and n±{(p)]. More quantitatively, one may compute the expectation value and the 
standard deviation for the number of atoms with a given spin projection in a condensate 
either by using quantum mechanics with the operators no{(p), etc., or by averaging over 
quantum trajectory simulations of a large number of condensates. We have found that, 
for large N, the results agree. 
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5. Concluding remarks 



Both quantum trajectory simulations and (|I2|) give seemingly sensible predictions for 
the frequencies of spin projections and their fluctuations even for the en masse Stern- 
Gerlach experiment. But they do a lot more, too. For instance, a histogram such 
as in figure |^ may be prepared for the results of en masse experiments and compared 
with the predictions of quantum trajectory simulations. Such a histogram probes the 
spin statistics in principle to all orders, beyond average and standard deviation. We 
thus have a nontrivial prediction from ([T^ ) that may be contrasted with experiments. 
Another possible extension of our method is time evolution, as in the Josephson effect [^ 
or in the time dependence of spin correlations. We may straightforwardly incorporate 
time ordering and time evolution of the annihilation operators into (0). In contrast, an 
approach in terms of the number operators no{t), n±{t) without an explicit consideration 
of the measurements probably makes a subtle affair ||12|| . 

It would seem worthwhile to derive (|12D from microscopic arguments |15| , along the 



lines of earlier work on photons [|14|. To the extent that this can be done, our quantum 



trajectory simulations would be proven to be a universal method for the studies of the 
detected properties of a spinor condensate. More generally, coherence theory of light 
is built entirely on the theory of photon detection, so it appears plausible that the 
eventual coherence theory for Bose-Einstein condensates is based on detection theory as 
well. The combination of measurement theory and accompanying quantum trajectory 
simulations is at present our best bet for a general approach to the coherence properties 
of a condensate. 

In sum, we have laid down a measurement theory and the accompanying quantum 
trajectory simulations for a spinor condensate, and proposed a nontrivial experiment in 
which our approach may be tested. The main result thus far is that the state of the spinor 
condensate and measurements thereof are intertwined to the extent that it is probably 
impossible to distinguish experimentally between the symmetry-broken coherent state 
and the fragmented state of the condensate. In the process we have uncovered new 
aspects of measurement theory which should prove fruitful in future discussions of the 
coherence properties of a condensate. 
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